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^  ABSTRACT 

The  performance  of  aircraft  engines^ operating  in  areas  where  the 
atmosphere  is  polluted  by  small  solid  particles^can  suffer  due  to  blade 
surface  erosion.  ^The’^rosion  damage  can  lead  to  significant  reduction  in 
engine  efficiency  as  well  as  performance,  due  to  the  change  in  blade 
surfaces,  tip  leakages  and  blade  pressure  distribution.  This  report 
presents  the  results  of  an  investigation  of  the  solid  particle  dynamics  and 
the  resulting  blade  erosion  through  a  helicopter  engine  with  inlet  particle 
separator.  -The  particle  trajectories  are  computed  in  the  inlet  separator 
which  is  characterized  by  considerable  hub  and  tip  contouring  and  radial 
variation  in  the  swirling  vane  shape.  The  nonseparated  particle 
trajectories  are  determined  through  the  deswirling  vanes  and  the  five  stage 
axial  and  one  stage  radial  compressors.  ^The  ^impact  data  for  a  very  large 
number  of  ingested  particles  is  used  to  calculate  the  resulting  blade 
surface  erosion.  The  erosion  pattern  indicates  the  location  of  maximum 

A  / 

blade  erosion.  In  addition,  the  distribution  of  particle  impact  data/  which  ' 
-art  pertinent  to  erosion  such  a3  ‘he  impact  velocities,  impact  angles  and  in 
particular,  the  frequency  of  particle  impacts ,  ar^presented  to  suggest 
possible  procedures  to  reduce  the  erosion  in  the  critically  affected  blade 
areas .  ■"  '  ■  .'1  :  < „  c 
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NOMENCLATURE 

particle  drag  coefficient 
particle  diameter 

force  of  interaction  between  the  gas  and  the  particle 
meridional  distance 

radial  distance  from  the  turbomachine  axis 

velocity  relative  to  the  blade  surface 
axial  coordinate 

particle  impingement  angle  relative  to  the  blade  surface 
frequency  of  particle  impacts  (total  number  of  particle  impacts 
per  unit  area  of  the  blade  surface  per  unit  number  of  ingested 
particle) 

blade  surface  erosion  parameter,  mg  of  eroded  material  per  unit 
area  of  the  blade  surface  per  unit  mass  of  ingested  particles 
2 

(mg/m  /gm) 
angular  coordinate 
gas  density 

particle  material  density 


rotor  speed  (radians/s) 


Subscripts 


refers  to  particle  impact  conditions  relative  to  the  surface 
particle 

component  in  radial  direction 

refers  to  particle  rebound  conditions  relative  to  the  surface 
component  in  axial  direction 
component  in  circumferential  direction 
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INTRODUCTION 


The  problem  of  a  gas-flow  mixed  with  solid  particles  in  axial  and 
radial  flow  turbomachines  has  great  importance  in  aeronautical,  industrial, 
and  naval  applications.  The  performance  of  the  aircraft  engines  is  known  to 
deteriorate  rapidly  when  they  operate  in  areas  where  the  atmosphere  is  laden 
with  solid  particles.  In  general,  these  particles  do  not  follow  the  flow 
streamlines  in  the  engine  and  tend  to  impact  the  blade  surfaces  due  to  their 
higher  inertia.  These  blade  impacts  can  cause  severe  erosion  damage  that  i3 
manifested  in  the  pitting  and  cutting  of  the  leading  and  trailing  edges,  and 
a  general  increase  in  the  blade  surface  roughness  [1],  Excessive  erosion 
damage  can  lead  to  the  structural  failure  of  the  blades. 

The  erosion  problems  in  military  and  commercial  airplane  gas  turbine 
engines  are  generally  recognized  in  industry.  The  operating  life  of 
helicopter  engines  performing  in  sandy  areas  is  very  short,  from  50  to  250 
hours  [2].  Recently,  many  reported  industrial  gas  turbine  engine  failures 
have  been  connected  to  their  exposition  to  particulate  flows  [3].  Research 
efforts  are  still  needed  to  provide  a  thorough  knowledge  of  the  various 
parameters  which  influence  the  extent  of  erosion  damage. 

The  erosion  of  metals  by  solid  particle  laden  flow  have  been 
investigated  experimentally  [4,  5,  6]  to  determine  the  parameters  that 
affect  the  erosion  rate  and  to  arrive  at  empirical  equations  for  the  target 
material  erosion  in  terms  of  these  parameters.  The  results  or  these 
investigations  demonstrate  that  for  a  given  particle-target  material 
combination,  the  erosion  rate  is  affected  by  the  impacting  velocity, 
impingement  angle,  and  by  the  metal  and  the  gas  temperatures.  Prediction  of 
turboraachine  blade  erosion  hence  requires  that  a  large  number  of  particle 
trajectory  calculations  are  carried  throughout  the  machine  to  determine  the 
pattern  and  frequency  of  particle  blade  impacts  as  well  as  the  magnitude  and 
direction  of  the  impacting  velocities  relative  to  the  blade.  This  data  can 
then  be  combined  with  the  empirical  equation  to  calculate  the  blade  erosion 
pattern  and  intensity. 

Erosion  cascade  tunnels  have  also  been  used  to  determine  experimentally 
the  change  in  compressor  cascade  performance  due  to  erosion  [7].  The 
changes  in  the  airfoil  configuration  and  its  surface  roughness  were  found  to 
influence  both  the  blade  loading  and  the  cascade  1033  coefficient.  However, 
blade  erosion  in  multistage  turbomachines  differ  significantly  from  the 
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cascade  erosion  due  to  the  difference  in  the  particle  impact  locations  as 
well  as  in  their  impacting  velocities  and  impingement  angles  relative  to 
these  surfaces.  Since  the  particles  impact  the  blades  as  well  as  the 
casings,  the  geometrical  configuration  of  the  latter  would  also  influence 
the  blade  erosion.  It  was  shown  that  the  hub  contour  geometry  [8]  and  the 
radial  variation  in  the  blade  shape  [9]  not  only  affect  the  particle 
trajectory  but  also  their  radial  and  circumferential  distribution  after  the 
blade  row.  Furthermore ,  since  the  particles  are  generally  redistributed  in 
the  radial  direction  under  the  influence  of  the  centrifugal  forces  as  they 
travel  through  the  machine  and  repeatedly  impact  the  different  blade 
rows  [10,  11,  12],  the  erosion  of  the  blades  in  the  multistage  machines  is 
dependent  on  their  stage  locations. 

This  work  presents  the  results  of  a  detailed  study  of  blade  and  vane 
surface  erosion  through  the  inlet  separator,  in  the  deswirling  vanes,  and 
the  five  stage  axial  compressor,  and  the  radial  compressor  of  the  G.E.  T700 
helicopter  gas  turbine  engine  described  in  reference  [13].  The  inlet  is 
characterized  by  swirling  vanes  placed  in  a  highly  contoured  hub  zone  to 
centrifuge  the  particles  for  separation.  The  unscavanged  flow  proceeds 
through  de3wirling  vanes  placed  in  an  annulus  zone  with  decreasing  inner  and 
outer  radii  before  the  five  stage  axial  flow  compressor. 

Two  separate  computer  codes,  one  for  predicting  the  particle  dynamics, 
and  another  for  predicting  the  erosion  are  used  to  obtain  the  presented 
results.  The  three  dimensional  particle  trajectory  calculation,  in  the 
various  components  of  the  machine  [12,  14,  15,  16],  provide  the  particle 
blade  impact  locations  and  the  magnitude  and  direction  of  the  particle 
impacting  velocity  relative  to  each  blade  [3].  A  computer  code  that  is 
especially  developed  to  predict  the  blade  surface  erosion  uses  the  computed 
blade  impact  data  and  the  experimental  results  for  blade  material  erosion. 
The  presented  results  show  the  magnitude  and  pattern  cf  blade  erosion  in  the 
various  blade  row3.  The  results  indicate  that  the  deswirling  and  swirling 
vane  erosion  is  strongly  influenced  by  the  inlet  separator  geometry.  The 
multistage  axial  compressor  blade  erosion  is  found  to  be  strongly  influenced 
by  the  blade  row  location.  The  results  also  indicate  that  the  erosion 
intensity  in  the  compressor  increases  in  the  latter  stages  and  that  higher 
local  values  of  erosion  are  generally  encountered  in  the  rotor  Hades  than 
in  the  stator  blades.  The  results  of  the  radial  compressor  indicate  that 
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the  impeller  blade  pressure  surface  encounters  on  erosion  intensity  higher 
than  that  of  the  splitter  surfaces.  It  was  also  found  that  the  erosion 
intensity  of  the  radial  impeller  is  as  high  as  the  axial  compressor  later 
stages. 


ANALYSIS 


The  equations  governing  the  particle  notion,  in  the  turboaachinery  flow 
field  are  written  in  cylindrical  polar  coordinates  relative  to  a  frame  of 
reference  fixed  with  respect  to  the  rotating  blade. 


d0n 

F  ♦  (—£• 
r  vdt 


(1) 


de 

(__JE  ♦  w) 
dx 


(2) 


where  r  ,  ,  z^  define  the  particle  location  in  cylindrical  polar 

coordinates,  and  c»  the  blade  angular  velocity.  The  centrifugal  acceleration 
and  Corioli3  acceleration  are  represented  by  the  last  terms  on  the  right 
hand  side  of  equations  (1)  and  (2),  The  first  term  on  the  right  hand  side 
of  equations  (1)  through  (3)  represents  the  force  of  interaction  between  the 
two  phases,  per  unit  mass  of  particles.  This  force  is  dependent  on  the 
relative  velocity  between  the  particles  and  the  gas  flow,  as  well  as  the 
particle  size  and  shape.  Under  the  particulate  flow  conditions  in 
turbomachines,  the  effect  of  the  forces  due  to  gravity  and  to  interparticle 
interactions  are  negligible  compared  to  those  due  to  the  aerodynamic  and 
centrifugal  forces.  In  addition,  the  force  of  interaction  between  the  two 
phases  i3  dominated  by  the  drag  due  to  the  difference  in  velocity  between 
the  solid  particles  and  the  flow  field.  The  force  of  interaction  per  unit 
mass  of  solid  particles  is  given  by 
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where  V  ,  V„,  V  ,  represent  the  relative  gas  velocities  in  the  radial, 
r  v  z 

circumferential  and  axial  directions,  respectively  and  p,  p  are  the  ga3  and 


^  solid  particle  material  densities,  d  the  particle  diameter,  and  CQ  the 

,  particle  drag  coefficient.  This  coefficient  is  dependent  on  the  Reynolds 

ft  number,  which  is  based  on  the  relative  velocity  between  the  particle  and  the 

ga3.  Empirical  relations,  as  shown  in  reference  [14],  are  used  to  fit  the 
y,  drag  curve  over  a  wide  range  of  Reynolds  numbers. 


% 


£ 


■r: 


I 


V 


& 


.  / 
.V 
.V 


I  - 


Trajectory  Calculations 

The  particle  trajectory  calculations  consist  of  the  numerical 
integration  of  the  equations  (1)-(4)  in  the  flow  field,  up  to  the  point  of 
blade,  hub  or  tip  impact.  The  magnitude  and  direction  of  particle 
rebounding  velocity  after  these  impact?  are  dependent  on  the  impacting 
conditions  and  the  particular  particle  and  surface  material  combination 
under  consideration. 

The  particle  rebounding  conditions  are  determined  from  empirical 
correlations  of  the  restitution  parameters.  These  correlations  [17]  are 
based  on  the  experimental  data  obtained  using  LDV  for  particle  laden  flows 
over  metal  samples  at  various  incidence  angles  and  flow  velocities  in  a 
special  tunnel  [5,  6].  The  ratios  were  found  to  be  mainly  dependent  upon 
the  impingement  angle  [17]  for  a  given  particle-material  combination. 

In  the  present  3tudy  the  magnitude  and  direction  of  the  particle 
rebounding  velocity  after  impact  were  obtained  using  the  following  empirical 
equation  which  was  obtained  from  [17]  for  Kingston  fly  ash  particles 
impacting  on  a  410  stainless  steel  target.: 

VP  /VP  “  1  ’  0-90847  3t  +  0.3072  Bj  *  0.05695  8j 
R  I 

(5) 


and 


eD/eT  «  i  -  0.38746  eT  +  0.51442  b^  ♦  0.45094  e? 

n  1  1  1  l 


(6) 


where  Vp  and  8.  represent  the  particle  impact  velocity  and  impingement 
PI  1 


angle,  while  Vn  and  Ba  are  the  rebound  velocity  and  rebound  angle  measured 

K 


from  the  impacted  surface. 

In  the  case  of  rotor  blade  impacts,  the  velocities  and  angles  of  impact 
and  rebound  are  considered  relative  to  the  rotating  blade.  Furthermore,  in 
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order  to  use  the  empirical  correlations  obtained  from  two  dimensional 
measurements,  in  the  three  dimensional  particle  trajectories,  a  knowledge  of 
the  geometry  of  the  solid  surface  at  the  impact  location  is  required. 
Geometrical  data  describing  the  blade,  hub  and  tip  configuration  need  to  be 
available  during  the  trajectory  calculations  for  use  in  the  determination  of 
the  impact  location  and  the  impingement  angle  relative  to  the  surface. 

Blade  Passage  Geometric  Representation 

The  accurate  representation  of  the  blade  passage  geometry  is  civ^ial  in 
the  particle  trajectory  calculations  and  blade  erosion  prediction  as  it 
influences  the  predicted  blade  impact  locations  and  particle  impact  velocity 
and  impingement  angle  relative  to  the  blades.  In  addition,  the  particle 
trajectory  after  each  impact  is  also  affected  by  these  impact  conditions 
through  the  rebound  restitution  correlations  of  equations  (5)  and  (6).  Even 
if  the  particle  distribution  is  uniform  as  particles  enter  the  blade  row, 
they  are  generally  redistributed  in  a  nonuniform  pattern  as  they  leave  due 
to  blade  surface  impacts  [16],  The  blade  airfoil  shape  and  inner  and  outer 
casing  geometry  must  therefore  be  accurately  represented  in  the  trajectory 
analysis  as  it  strongly  influences  the  predicted  particle  dynamics.  In 
addition,  the  radial  variation  in  the  blade  shape  Is  also  important  as  it 
affects  the  particle  rebound  characteristics  in  the  three  dimensional 
trajectories  [33.  and  might  influence  the  particle  redistribution  in  the 
radial  direction.  The  inner  and  outer  casing  coordinates  and  the  blade 
surface  coordinates  at  several  sections  between  the  hub  and  tip  are  used  in 
the  trajectory  calculations  to  determine  the  blade  surface  coordinates  at 
other  locations  using  interpolation.  At  least  two  blade  sections  are 
required  but  more  sections  can  be  used  to  describe  the  hub  to  tip  variation 
in  the  blade  3hape. 

Flow  Field  Representation 

The  flow  field  affects  the  particle  trajectories,  through  the 
aerodynamic  forces  of  interaction  consisting  mainly  of  the  drag  force. 
Pitchwise  as  well  as  radial  and  axial  flow  field  variations  must  be 
represented  in  the  particle  trajectory  calculation.  Careful  consideration 
must  be  given  to  the  flow  field  description  in  the  particle  trajectory 
calculations  [8],  in  order  to  represent  the  important  flow  field 
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characteristics  without  excessive  computer  storage  requirements.  The  radial 
variation  in  the  flow  field  is  affected  by  the  radial  variation  in  blade 
loading  and  by  inner  and  outer  casing  radii  contours. 

Referring  to  Fig.  1 ,  one  can  observe  that  the  inner  and  outer  casings 
are  highly  contoured  in  the  inlet  separator  and  through  the  deswirling 
vane3,  than  through  the  five  stage  axial  flow  compressor.  The  flow  field 
representation  in  the  particle  trajectory  calculation  was  hence  chosen 
accordingly  to  reflect  the  important  three  dimensional  flow  characteristics. 
The  flow  field  on  a  mean  hub  to  tip  stream  surface  and  a  blade  to  blade 
velocity  gradient  Cl 8]  are  used  to  represent  the  flow  field  in  the  particle 
trajectory  code  for  the  inlet  separator  with  the  swirling  vanes  [16]  and  in 
the  deswirling  vanes.  On  the  other  hand,  in  the  axial  compressor  particle 
trajectory  calculations,  the  flow  field  in  the  compressor  blade  passage  is 
synthesized  from  a  number  of  blade-to-blade  stream  surface  solutions  Each 
flow  field  solution  on  a  blade-to-bxade  stream  surface  [18]  provides 
accurate  description  of  the  pitchwise  and  streamwise  variation  in  the  flow 
properties.  The  stream  surface  shapes  and  stream  filament  thickness  as 
determined  from  a  mid  channel  hub  to  tip  stream  surface  solution  [19]  model 
the  influence  of  the  inner  and  outer  casing  contours  and  the  hub  to  tip 
variation  in  the  blade  shapes. 

Calculation  of  the  three  dimensional  flow  field  through  the  radial 
compressor  impeller  was  conducted  in  two  steps  [20].  The  first  3tep  uses 
the  code  of  Ref.  [18]  to  get  the  flow  in  the  hub  shroud  meridional  plane  and 
the  stream  tube  filament  thickness  which  are  the  input  parameter  for  the 
second  step.  The  second  step  uses  the  panel  method  [21]  to  obtain  the  blade 
to  blade  solutions  on  a  number  of  stream  surfaces  stacked  from  hub  to 
shroud. 

Particle  Trajectory  Computation  Procedure: 

One  can  conclude  from  the  previous  discussion  that  the  spanwise, 
pitchwise  and  chordwise  variations  in  flow  properties  are  therefore 
represented  in  the  particle  trajectory  calculations.  In  addition  the  effect 
of  the  radial  variation  in  the  blade  shape  and  hub  and  tip  contouring  is 
taken  into  consideration  in  the  particle  trajectory  calculations.  The 
output  of  the  particle  trajectory  calculation  code  consists  of  three  files. 
The  first  file  provides  the  input  for  particle  trajectory  plots.  The  second 


particle  trajectory  output  file  contains  all  particle  exit  conditions  for  a 
given  blade  row  and  is  used  as  input  to  provide  the  initial  conditions  for 
the  particle  trajectories  in  the  following  blade  row.  This  data  is  first 
adapted  to  the  new  frame  of  reference,  which  is  fixed  in  each  blade  row. 

The  particle  axial  and  radial  location  is  not  affected  by  the  relative 
motion  between  successive  blade  rows.  The  circumferential  variation  in  the 
particle  distibution  at  the  exit  of  a  given  blade  row,  is  in  reality  evened 
out  at  the  inlet  to  the  following  blade  row  due  to  the  relative  motion 
between  successive  blade  rows.  This  effect  is  simulated  in  the  computations 
by  re-distributing  the  particles  randomly  in  the  circumferential  direction, 
in  the  new  frame  of  reference  which  is  fixed  relative  to  the  subsequent 
blade  row.  The  third  file  contains  all  computed  particle  impact  data, 
including  impact  locations,  and  the  impact  velocity  and  impingement  angle 
relative  to  the  solid  surfaces,  for  use  in  the  code  for  blade  surface 
erosion  computations. 

Blade  Mass  Erosion  Prediction 

The  computation  of  the  blade  surface  erosion  combines  the  empirical 
correlations  for  the  mass  erosion  parameter  for  a  given  particle  target 
material  combination  [4,  5]  with  the  results  of  the  particle  trajectory 
calculations  for  a  large  number  of  particle  trajectories  in  the  form  of  the 
impact  velocities  and  impingement  angles. 

The  experimental  measurements  in  the  erosion  tunnel  indicate  that  the 
target  erosion  was  found  to  be  dependent  upon  the  particle  impact  velocity, 
impingement  angle  and  on  the  flow  and  target  temperatures.  Experimental 
measurements  were  obtained  in  the  erosion  tunnel  for  quartz  particles  (165 
micron  diameter)  impacting  510  stainless  steel  at  different  impacting 
velocities  and  impingement  angles. 

The  experimental  data  was  used  to  obtain  the  following  empirical 
equation  for  the  erosion  mass  parameter,  e,  which  is  defined  as  the  ratio  of 
the  eroded  mass  of  the  target  material  to  the  mass  of  impinging  particles. 

QO  2  2  2  2 

e  -  K1  [1  ♦  CK(K12  sin(j-  B,))]  ^  ]  cos  ^  (1  -  Rt) 

o 

♦  K3  (V1  sing,)4 

where  and  B^  are  the  impact  velocity  and  impingement  angle,  respectively. 


Rt  -  1  -  0.0016  Vj  »In81 

Bq  -  angle  of  maximum  erosion 

CK  -  1  for  6,  <  2B 
i  —  o 

CK  -  0  for  81  >  2Bq 

K1  ,  ^  and  are  mater'ial  constants  which  were  found  to  be  as  follows: 

K1  -  0.5225  x  10-5 

K12  -  0.266799 

K3  «  0.5^9  x  10~12 

The  erosion  computer  code  was  especially  developed  to  calculate  the 
blade  surface  erosion  distribution.  The  program  input  consists  of  the  blade 
surface  geometrical  data  and  the  particle  blade  impact  data  as  determined  by 
the  particle  trajectory  calculations.  A  blade  surface  fitted  grid  is 
generated  and  the  impact  data  is  then  combined  with  the  empirical  equations 
to  predict  the  mass  erosion  intensity  and  pattern.  The  mesh  size  was 
carefully  selected,  according  to  the  finite  number  of  particles 
trajectories,  to  give  the  most  accurate  representation  of  the  blade  surface 
erosion. 

The  output  from  this  program  Includes  the  blade  erosion  distribution  as 
well  as  the  distribution  of  the  impact  velocity,  impingement  angles  and  the 
frequency  of  particle  impacts  on  the  blade  surface.  The  detailed 
distribution  of  all  these  significant  erosion  parameters  over  each  blade 
surface  in  a  multistage  compressor  is  essential  for  taking  specific  measures 
to  alleviate  blade  erosion  problems. 


RESULTS  AMD  DISCUSSION 

The  results  of  the  particle  trajectory  computations  and  blade  surface 
erosion  analysis  are  presented  for  165  micron  diameter  quartz  particles 
through  the  inlet  separator,  the  axial  and  the  radial  flow  compressors  of 
the  G.E.  T700  helicopter  engine  which  is  shown  in  Fig.  1.  Three  sets  of 
results  are  presented  in  each  component,  the  first  shows  the  particle 
trajectory  projections  on  the  blade  to  blade  surface  and  on  a  hub  to  tip 
surface.  The  second  set  of  data  is  also  pertinent  to  the  particle 
trajectory  computations  and  shows  the  actual  locations  of  particle  blade 
surface  impacts.  The  third  set  presents  the  results  of  blade  surface 
erosion  analysis  in  the  form  of  the  distribution  over  the  blade  surface,  of 
the  frequency  of  particle  impacts,  the  magnitudes  of  the  impacting  velocity 
and  impact  angle  relative  to  the  blade,  and  finally  the  predicted  blade 
erosion. 

Inlet  Separator 

Figure  1  shows  the  geometry  of  the  inlet  separator  which  combines  a 
highly  contoured  inner  casing  with  the  swirling  vanes  to  deflect  the 
particles  in  the  radial  direction.  The  particle  initial  conditions  were 
taken  statistically  to  represent  uniform  particle  loading  with  uniform  axial 
velocity  upstream  of  the  engine  inlet.  The  particle  trajectory  calculations 
were  performed  in  the  inlet  separator  as  described  in  reference  [16],  The 
results  of  the  particle  trajectory  computations  are  presented  in  figures  2 
and  3  which  show  the  projection  of  sample  particle  trajectories  in  the  blade 
to  blade  cylindrical  surface  and  the  meridional  plane  through  the  machine 
axi3,  respectively.  One  can  observe  that  the  trajectories  of  these  large 
particles  deviate  considerably  from  the  flow  streamlines  because  of  their 
high  inertia.  These  particles  can  be  seen  to  continue  in  their  axial  motion 
uninfluenced  by  the  flow  field  until  they  impact  the  blade  or  casing.  One 
can  see  that  some  of  the  particles  near  the  tip  traverse  the  twisted  blades 
without  any  impacts.  Many  of  these  particles  are  still  separated  however 
because  of  their  outward  radial  location  relative  to  the  splitter  between 
the  engine  and  scavanged  flow.  The  particles  that  impact  the  swirling  vanes 
pressure  surface,  rebound  with  a  smaller  radially  outward  velocity  component 
and  a  larger  circumferential  velocity  component  that  causes  their 
centrifugion.  These  particles  impact  the  tip  casing  and  acquire  a  radially 
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inward  velocity  after  rebounding  which  causes  thee  to  enter  the  engine.  On 

the  other  hand,  many  of  the  particles  that  impact  the  inner  casing  are 
reflected  in  the  radially  outward  direction  towards  the  tip  annulus.  Many 
of  these  particles  also  impact  the  outer  casing  and  acquire  a  radial  inwarad 
velocity  after  rebounding  and  later  enter  the  engine.  The  very  few 
particles  exempt  from  this  behavior  are  those  impacting  the  inner  casing 
further  downstream,  much  closer  to  the  hump,  and  might  miss  entering  the 
engine  due  to  their  later  impact  locations  with  the  outer  casing,  which 
causes  them  to  separate. 

Figures  4  and  5  show  the  particle  impact  locations  on  the  swirling 
vanes  suction  and  pressure  surfaces  respectively.  It  is  obvious  from  these 
figures  that  the  particles  impact  the  vane  pressure  surface  much  more 
frequently  than  the  suction  surface.  Figure  5  also  shows  that  the  particle 
impacts  with  the  vane  increases  towards  the  hub  and  the  trailing  edge. 

The  particle  impact  frequency  impinging  velocity  and  angle  distribution 
over  the  vane  surface  was  determined  from  a  large  number  of  particle 
trajectories  and  the  results  are  3hown  in  Figs.  6  through  8  respectively  for 
the  pressure  surface.  Figure  6  confirms  that  generally  the  particle  impact 
frequency  increases  towards  the  hub  except  at  the  trailing  edge.  A  high 
intensity  impact  region  extends  over  most  of  the  vane  chord  towards  the  hub. 
Two  other  regions  of  high  impact  frequency  are  also  predicted  at  the  corner 
of  the  trailing  edge  and  outer  casing,  and  the  opposite  corner  of  the 
leading  edge  and  Inner  casing.  Figure  7  shows  the  distribution  of  the 
particle  impacting  velocity  on  the  pressure  surface  of  the  swirling  vanes. 
The  highest  impacting  velocities  (60  m/s  or  more)  can  be  seen  in  a 
triangular  part  of  the  blade  near  the  leading  edge  and  the  tip  which  extends 
up  to  75t  of  both  the  span  and  chord.  This  corresponds  to  the  particle 
inlet  velocities,  which  remain  practically  unchanged  until  the  first 
particle  impact.  The  low  impact  velocities  in  the  other  parts  of  the  blade 
can  be  attributed  to  the  lower  velocity  of  many  of  the  particles,  impacting 
thi3  area  after  encountering  previous  impacts  with  the  blade  surface,  the 
hub  or  tip  casing  which  causes  a  loss  in  their  momentum.  The  higher 
impacting  velocity  regions  near  the  trailing  edge  reflect  first  time  impacts 
for  the  particles  that  traverse  roost  of  the  blade  chord  with  no  previous 
interactions  to  their  blade  impacts  near  the  trailing  edge.  Figure  8  shows 
the  impact  angle  distribution  over  the  swirling  vanes  pressure  surface.  The 
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highest  impact  angles  are  seen  at  the  leading  edge  area  which  then  decreases 
towards  the  mid  chord,  and  then  increases  again  towards  the  trailing  edge  as 
the  blade  angle  increases. 

An  erosion  code  was  especially  developed  to  predict  the  blade  surface 
erosion  based  on  the  blade  surface  impact  data  for  a  large  number  of 
particles  and  the  empirical  equations  from  the  experimental  measurements  for 
the  particle  target  material.  Figure  9  shows  the  distribution  of  the 
calculated  erosion  parameter  on  the  pressure  surface  as  determined  from  a 
large  number  of  particle  trajectories.  This  figure  shows  that  the  maximum 
blade  erosion  is  predicted  at  the  corner  of  the  blade  between  the  trailing 
edge  and  the  hub.  The  particle  impact  angle  of  25-30  degrees  in  this  region 
is  close  to  the  impact  angle  corresponding  to  the  maximum  erosion  for 

this  particle  target  material  combination.  Both  impact  velocity  and 
frequency  of  particles  impact  are  also  high  in  this  region. 


Deswirling  Vane3 


The  unseparated  particles  enter  the  engine  with  the  unscavenged  air 
which  follows  the  inner  part  of  the  flow  path  through  the  separator. 

Figures  10  and  11  present  sample  particle  trajectories  through  the 
deswirling  vanes  viewed,  in  the  z-0  plane  and  in  the  z-r  plane  respectively. 
Many  of  the  unseparated  particles  enter  with  very  small  radial  velocities 
and  consequently  impact  the  contoured  outer  casing.  Some  of  these  and  the 
other  particles  impact  either  the  pressure  side  or  the  suction  side  of  the 
vanes. 

Figures  12,  13  show  the  locations  of  particle  impacts  with  the 
deswirling  vane  surfaces.  One  can  3ee  that,  over  the  vane  suction  surface, 
the  particle  impacts  are  confined  to  one  third  the  vane  chord  towards  the 
leading  edge,  and  in  another  narrower  band  that  extends  from  the  leading 
edge  blade  corner  at  the  tip  to  about  one  quarter  the  blade  span  at  the 
trailing  edge.  On  the  other  hand  the  particle  impacts  with  the  pressure 
surface  are  concentrated  towards  the  outer  half  of  the  blade  span,  and 
extend  more  over  the  vanes  leading  edge  and  trailing  edge  and  less  in 
between.  These  particle  impact  patterns  are  a  result  of  the  combined  effect 
of  the  particle  distribution  as  they  enter  the  vanes  and  the  geometry  of  the 
twisted  vanes  and  the  highly  contoured  annulus. 
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Figures  1*)  through  17  present  the  distribution  of  the  erosion 
parameters  e,  V,  8  and  e  over  the  suction  3ide  of  the  deswirling  vanes.  One 
can  see  from  figure  1  *t  that  the  highest  values  of  particle  impact 
frequencies  are  encountered  near  the  tip  and  leading  edge  corner  and 
decrease  gradually  towards  the  hub  and  the  trailing  edge.  In  addition.  Fig. 
15  shows  that  the  impacting  velocities  are  highest  in  the  same  areas  of 
higher  impact  frequencies.  Figure  16  shows  that  the  particle  impact  angles 
are  generally  low  with  the  highest  values  of  impacting  angles  close  to 
&o“25°,  at  the  leading  edge  close  to  the  tip.  This  corresponds  to  the 

impact  angle  of  maximum  erosion.  All  three  effects  combined  together  result 
in  the  highest  erosion  rate  on  the  blade  suction  surface  at  the  leading  edge 
tip  corner  as  shown  in  Fig.  17. 

Figures  18  through  21  present  the  erosion  parameters  distribution  on 
the  deswirling  vanes  pressure  surface.  Comparing  Figs.  18  and  1*1,  one  can 
see  that,  more  particles  impact  the  pressure  surface  than  the  suction 
surface.  The  particle  impact  frequency  distribution  is  al3o  different,  with 
maximum  values  on  the  pressure  surface  near  the  vane  mid  chord  at  the  tip. 
Figure  19  shows  that  the  particle  impact  velocities  do  not  change 
drastically  over  the  deswirling  vanes  pressure  side.  Comparing  Figures  19 
and  15,  one  can  see  that  the  particle  impacting  velocities  are  generally 
lower  over  the  pressure  surface  than  the  suction  surface.  Figure  20  shows 
that  the  particle  impact  angles  are  generally  higher  than  8o,  the  impact 

angle  of  maximum  erosion.  The  lowest  values  of  the  particle  impact  angles 
are  at  the  vane  tip  near  the  leading  corner  where  it  i3  close  to  the  value 
of,  Bo»25°,  the  angle  of  maximum  erosion.  Figure  21  shows  that  the  erosion 

pattern  on  the  pressure  side  of  the  deswirling  vanes  is  somewhat  similar  to 
that  of  the  impacting  frequency,  with  the  maximum  at  the  tip  of  the  blade 
near  the  mid-chord.  Comparing  Figs.  21  and  17  one  can  see  that  the  maximum 
value  of  the  deswirling  vane  pressure  surface  erosion  is  almost  the  same  as 
for  the  suction  surface,  in  spite  of  the  much  higher  impacting  frequencies 
over  the  pressure  surface.  Thi3  can  be  attributed  to  the  higher  impacting 
velocities  on  the  suction  surface,  Fig.  15. 


Axial  Compressor 


The  blade  erosion  characteristics  of  the  G.E.  T700  helicopter  engine 
five  stage  axial  flow  compressor  were  studied  on  the  basis  of  the  trajectory 
data  presented  in  reference  [33.  The  particle  trajectory  calculation,  for 
165  microns  diameter  particles,  were  performed  throughout  the  five  stages 
axial  flow  compressor.  The  results  of  the  first  and  last  stages  particle 
trajectory  are  shown  in  Figs.  22  through  29.  respect ively.  One  can  notice 
from  Figs.  23,  25  and  27  that  only  the  pressure  side  of  the  rotor  and  stator 
blades  are  subjected  to  particle  impacts.  Details  of  the  trajectory  results 
and  the  location  of  particle  impacts  on  the  blade  surfaces  can  be  found  in 
reference  [33. 

The  data  required  by  the  code  which  was  specially  developed  for 
analyzing  the  erosion  characteristics  on  the  machine  blade  surfaces 
consisting  of  the  impact  characteristics  of  a  large  number  of  particles  wa3 
obtained  from  the  trajectory  computations.  Only  the  blade  pressure  surfaces 
data  was  considered  because  the  blade  suction  surface  erosion  is 
insignificant.  Figures  30  through  3^  present  the  important  blade  erosion 
parameters,  namely,  particle  impacting  velocity  V,  particle  impingement 
angle  6,  the  frequency  of  particle  impacts  9,  and  particle  erosion  parameter 
e,  for  the  five  rotors,  of  the  axial  compressor.  One  can  observe  from  these 
figures  that  both  particle  impact  frequency  and  impact  velocity  patterns 
over  the  rotor  pressure  surface  are  somewhat  similar.  The  iapingeaent  angle 
pattern  is  different,  and  changes  from  a  maximum  value  at  the  hub  to  a 
minimum  near  the  tip  towards  the  trailing  edge.  The  blade  erosion  pattern 
i3  similar  to  both  particle  impact  frequency  and  impact  velocity  patterns 
with  the  maximum  erosion  at  the  blade  leading  edge  tip  corner,  then 
diminishing  towards  the  trailing  edge  hub  corner.  A  region  of  the  rotor 
blade  pressure  surface  at  the  hub  trailing  edge  corner  extending  diagonally 
over  one  third  of  the  blade  chord  near  the  hub  to  about  two  thirds  of  the 
blade  span  at  the  trailing  edge  is  not  subjected  to  any  particle  impacts  and 
hence  suffer  no  erosion. 

Figures  35  through  39  show  the  variation  of  the  erosion  parameters  v, 

8,  9,  and  e  on  the  pressure  surfaces  of  the  five  stators  of  the  axial  flow 
compressor.  One  can  see  from  these  figures  that  the  stator  blade  pressure 
surface  erosion  pattern  i3  very  similar  to  the  particle  impact  frequency 
pattern.  The  maximum  stator  blade  erosion  is  always  at  the  tip  but  is  not 


necessarily  near  the  leading  edge,  A  small  portion  of  the  stator  blade  at 
the  leading  edge  near  the  hub  is  not  subjected  to  particle  impacts.  This  is 
due  to  the  effect  of  the  preceding  rotor  impacts  which  centrifuge  the 
particles  in  the  radially  outward  direction.  However,  the  particles  are 
distributed  more  uniformly  in  the  radial  direction  as  they  leave  the  stator 
blades.  This  explains  the  difference  in  the  location  of  the  non  impacted 
blade  surfaces  in  stationary  and  rotating  blade  rows.  In  general  the 
erosion  intensity  increases  in  the  latter  stages  due  to  the  reduced  flow 
path.  The  maximum  value  of  stator  blade  erosion  are  generally  lower  than 
the  maximum  rotor  blade  erosion. 
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Centrifugal  Compressor 

The  particle  trajectories  as  well  as  the  3tudy  of  the  erosion 
parameters  were  conducted  through  the  centrifugal  compressor,  in  order  to 
determine  the  particle  inlet  conditions  at  the  centrifugal  compressor.  The 
trajectory  calculations  for  the  165  micron  particles  were  continued  through 
the  annular  gap,  between  the  axial  and  the  centrifugal  compressors.  The 
particle  trajectory  computations  were  carried  out  in  thi3  ax i symmetric  flow 
field  as  determined  using  the  computer  code  of  reference  [18].  The 
centrifugal  compressor  impeller  configuration,  as  shown  in  Fig.  ** ,  has  two 
different  size  splitters  between  the  blades. 

Figure  Ala  and  Alb  show  a  sample  of  the  particle  trajectories  projected 
in  the  M-e  and  r-z  planes  respectively  through  the  impeller  of  the 
compressor.  One  can  also  see  that  the  particles  initially  have  a  small 
radial  inward  velocity  component  which  they  had  a3  they  left  the  fifth 
stator  of  the  axial  compressor.  Shortly,  the  particles  outward  radial 
migration  starts  and  is  more  pronounced  after  they  impact  either  the  blade 
or  the  hub  surface.  It  is  clear  from  Fig.  4la  that  the  particles  impact  the 
pressure  surfaces  of  the  impeller  blade  and  the  second  splitter.  Figure  41b 
shows  many  particle  impacts  with  the  hub  and  very  few  particle  impacts  ;with 
the  shroud. 

The  particle  impact  locations  with  the  blade  and  splitter  surfaces  are 
shown  in  Fig.  42.  One  can  see  from  this  figure  that  the  first  splitter 
suction  surface  is  subjected  to  some  particle  impacts,  while  most  of  the 
particle  impacts  are  with  the  second  splitter  pressure  surface  and  the  full 
blade  pressure  surface. 
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The  particle  impact  locations  for  the  impeller’s  hub  and  casing,  are 
presented  in  Fig.  ^3 -  Inspection  of  this  figure  shows  that  the  impeller  hub 
encounters  more  impacts  than  the  casing.  The  particle  distribution  at  the 
impeller  exit  is  given  in  Fig.  44,  which  shows  that  the  particles  are  more 
concentrated  near  the  pressure  side  of  the  impeller  blade  and  the  second 
splitter. 

The  conditions  of  impact  (velocity,  angle  and  location)  provided  by  the 
calculations  of  the  particle  trajectories  through  the  radial  impeller  were 
used  as  input.  An  erosion  study  of  the  radial  compressor  impeller  was 
conducted  using  the  same  procedures  outlined  previously  for  the  axial  flow 
compressor.  The  results  are  presented  for  the  erosion  parameters,  namely 
the  impact  velocity,  impact  angle,  impact  frequency  and  the  mass  erosion 
parameter  over  the  affected  blade  and  splitter  surfaces.  Figure  45  shows 
the  distribution  of  the  particle  impact  velocities  on  the  impeller  blade 
pressure  surface.  One  can  see  that  the  highest  impact  velocities  are  at 
inlet,  and  close  to  the  impeller  casing  at  outlet.  The  distribution  of  the 
impact  angles  on  the  same  blade  surface,  is  given  in  Fig.  46.  Inspection  of 
this  figure  shows  that  the  whole  blade  pressure  side  area  is  subjected  to 
impact  angles  that  correspond  to  maximum  erosion.  The  impeller  area  that  is 
subjected  to  the  highest  impacting  frequencies  constitutes  the  second  half 
of  the  impeller,  close  to  the  outlet  region,  as  shown  in  Fig.  47.  The 
resulting  blade  surface  erosion  over  the  blade  pressure  surface  is  given  in 
Fig.  48,  which  shows  high  values  of  erosion  at  the  impeller  inlet. 

The  distributions  of  the  erosion  parameters  over  the  pressure  surface 
of  the  second  splitter  vane  are  displayed  in  Figs.  49  through  52.  Figure  49 
shows  the  distribution  of  the  particle  impact  velocities  over  the  pressure 
side  of  the  second  splitter.  The  velocities  are  generally  high  comparable 
to  those  on  the  pressure  surface  of  the  impeller  blade.  The  highest 
velocities  are  close  to  the  impeller  casing.  The  distribution  of  the  impact 
angles,  Fig.  50,  show  that  most  of  the  surface  is  subjected  to  impact  angles 
close  to  that  corresponding  to  the  maximum  erosion.  The  areas  that  are 
subjected  to  the  maximum  frequency  of  impacts  are  the  same  areas  that  show 
the  maximum  rate  of  mass  erosion  as  shown  in  Figs.  51  and  52  respectively. 

Figures  53  through  56  show  the  distribution  of  the  same  erosion 
parameters  over  the  suction  side  of  tie  first  impeller,  where  the  particle 
impacts  are  only  limited  to  a  small  region  of  the  splitter  area  close  to  the 
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leading  edge,  (Fig.  3?).  Figure  53  shows  that  the  Impact  velocities  are 

relatively  low  compared  to  those  on  the  blade  pressure  surface  (Fig.  **5). 

The  impact  angle  distribution  over  the  same  splitter  surface  is  presented  in 

Fig.  5*4 ,  which  shows  that  smaller  impact  angles  than  that  corresponding  to 

maximum  erosion,  S  .  Figure  55  shows  the  area  of  the  splitter  vane  subject 
o 

to  the  highest  impact  frequencies.  Figure  56  shows  that  the  mass  erosion 
distribution  on  the  same  splitter  surface  is  similar  to  the  impact  frequency 
pattern  because  of  the  moderate  variation  in  both  impact  velocity  and  the 
impact  angle  distribution. 
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